Journal of Experimental Zoology Part B: Molecular and Developmental Evolution

High and low tides occur twice a day (every [~]12.4 hours), with the largest tidal ranges during spring tides around new and full moons (every [~]14.765 days). While these lunar cycles are known to influence many animal phenotypes, particularly the reproduction of coastal animals, the genetic basis of lunar-related rhythms remains unclear. Since phenotypic variation is a valuable resource for elucidating such mechanisms, we examined geographic variation in the lunar-regulated mass spawning of the grass puffer (Takifugu alboplumbeus) along the Japanese coast. We found that western populations spawn during the first half of the spring tides, whereas eastern populations spawn during the second half. Furthermore, although spawning typically occurs a few hours before high tide, this timing is restricted to a specific time window that is earlier in the western populations than in the eastern ones. Behavioral analysis of larvae also revealed a shorter free-running circadian period ({tau}) in the western population than in the eastern ones. As differences in {tau} affect individual variation in the timing of physiological functions and behaviors, we hypothesized that differences in {tau} could account for the different time windows and consequently the observed difference in spawning days. Population genomics analysis identified proline-rich transmembrane protein 1-like (prrt1l) as a candidate gene. Expression of prrt1l was observed in the circadian pacemaker suprachiasmatic nucleus, and triple CRISPR F0 knockout of prrt1l shortened the free-running period in larvae. These findings suggest a potential mechanism underlying the geographic variation in lunar-synchronized spawning behavior. HighlightsO_LIThe geographic variation exists in the lunar-regulated spawning of the grass puffer, with differences in spawning dates and times between western and eastern Japan. C_LIO_LIThe free-running period of western populations is shorter than that of eastern populations, which is consistent with their earlier spawning timing. C_LIO_LIPopulation genomics analysis identified prrt1l as a candidate gene harboring population-specific missense mutations, the knockout of which shortens the free-running period. C_LI

Dental enamel, the hardest mineralised tissue in the human body, has proven to be an excellent source of ancient proteins, which have been found to survive within dental enamel for at least twenty million years. In archaeological and palaeontological contexts, the enamel proteome is generally considered to be rather small, consisting of about twelve proteins, most of which are unique to enamel. During amelogenesis these proteins undergo in vivo digestion by matrix metalloproteinase 20 (MMP20) and kallikrein 4 (KLK4) as well as serine phosphorylation by family with sequence similarity member 20-C (FAM20C) that alter their characteristics. Gaining knowledge of the previously understudied influence of amelogenesis on the archaeological human dental enamel proteome could benefit various palaeoproteomic analysis, especially in an human evolutionary context. Here we present archaeological dental enamel proteomes and explore protein cleavage patterns and sequence coverage to estimate the effects of in vivo digestion, as well as explore phosphorylation patterns. Additionally, we present a new marker based on phosphorylation to estimate genetic sex.

In tunicates, larvaceans represent a fascinating case of evolution, where the chordate body plan has been maintained despite a rapidly evolving genome characterized by strong In contrast to other tunicates, larvaceans keep the chordate body plan during their entire life. They have acquired a highly specialized epithelium in charge of producing the "house", a complex extracellular apparatus used for filter feeding in the plankton. To what extent the house and this epithelium represent true molecular innovations withing chordates is a question for which thorough transcriptomics can bring novel insights. We conducted a developmental profiling of gene expression at the single-cell level in the larvacean Oikopleura dioica. We provide detailed descriptions of cellular transcriptomes associated with the house-synthesizing organ, which permits to define the molecular specifics of epithelial cell territories. We followed their emergence during development, and we identified genes that represent key candidate molecules for regulating the morphogenesis of the house-producing organ. Dynamic changes in gene expression and cell identities during major developmental transitions of the lifecycle illustrate that our dataset effectively allows access to the diversity of O. dioicas cell types in embryos and in adults. The resources presented here constitute critical assets to investigate larvacean biology and evolution for mechanistic and comparative goals.

Sarita Lake, British Columbia houses a distinctive population of threespine stickleback (Gastrosteus aculeatus L.) with a phenotype characterized by unusually large individuals relative to nearby conspecifics. We tested the hypothesis that members of this population are not isometrically larger but rather exhibit variation in allometric trajectories that reflect changes in developmental timing impacting the developmental-genetic architecture of the phenotype. We used 3D geometric morphometrics to characterize the size and shape of skulls, pectoral girdles and pelvic girdles from a sample of individuals from nearby freshwater and marine populations and compare them to a sample from Sarita Lake. We showed that individuals from the Sarita Lake population are larger in each body region compared to most other populations examined. Further, these individuals have dorsally expanded skulls and relatively robust pelvic armour. We also showed that the relationship between size and shape is differently structured among body regions and is heavily influenced by non-uniform sexually-mediated variation across populations sampled. Our results reflect complex underlying developmental trajectories, and we suggest that the large phenotype observed may be driven by fecundity selection on female size in combination with a limnetic trophic niche and relatively increased predation pressure in Sarita Lake.

Harbour porpoises (Phocoena phocoena) in the North and Baltic Seas are increasingly impacted by anthropogenic pressures, including underwater noise, fisheries and pollution. These pressures correlate with declining population health, particularly affecting the respiratory system. Growing pathological lesions, partly resulting from high prevalence of parasitic infestations and subsequent diseases, can impair tissue function and oxygen supply to distant end-organs. In this study, we applied an integrative MultiOmics approach (proteomics, metabolomics, lipidomics) to analyse the lungs and muscles of 12 wild harbour porpoises with compromised respiratory health. Our aim was to identify dysregulated biological pathways across omics layers to advance insights into adaptive physiological responses and to define disease-associated molecular signatures that could assist health assessments. Our analysis revealed pronounced immune system and antioxidative responses in the lungs and muscles, indicated by enhanced immunoglobulins, plasmalogens and glutathione-related proteins. In the lungs, high cardiolipin levels and reduced collagen suggest impaired tissue structure and function, while tissue maintenance processes were elevated in the muscle. Both tissues exhibited metabolic alterations suggestive of energetic imbalance, including increased purine metabolism in the lung and decreased lipid metabolism in the muscle. Several dysregulated molecules were shared across tissues, pointing to pathophysiological effects. The proposed disease-associated molecular signatures included the protein SLC25A4, the metabolite O-phosphoethanolamine and the lipid TG O-16:0_16:0_20:4 for the lung, and the protein SPEG, the metabolite pipecolic acid, and the lipid BMP 18:1_22:6 in the muscle. Our findings elucidate the complexity of molecular mechanisms linking anthropogenic and environmental stressors with vulnerability and resilience in a marine sentinel species. Furthermore, this study highlights the potential of integrative omics to define disease-related marker panels, thereby supporting ongoing and future health monitoring and conservation efforts.

The cockroach Blattella germanica possesses panoistic ovaries, in which oocytes lack nurse cells and therefore need to rely on their own transcriptional activity to support embryogenesis. Ovarian development in this species involves the development of a single basal ovarian follicle (BOF) per gonadotropic cycle, a process strictly regulated by endocrine signals, primarily juvenile hormone and ecdysone, which act at both the transcriptional and translational levels. In addition, transcriptional activity in these ovaries is necessary for both regulating and genome protection, and at this level, PIWI-interacting RNAs (piRNAs) play an essential role. Although insect ovaries are known to be particularly rich in piRNAs, their function in ovary maturation is still not well defined. For this purpose, we characterize the piRNA expression dynamics across seven key developmental and reproductive stages, ranging from late nymphal instars to post-vitellogenic adults. piRNA expression in B. germanica shows coordinated fluctuations. Expression remains stable in previtellogenic ovaries, whereas vitellogenic ovaries show pronounced changes. Moreover, vitellogenic ovaries exhibit reduced piRNA diversity due to strong enrichment of a subset of highly expressed piRNAs. Our data show that although piRNAs predominantly map to transposable elements, particularly LINEs, there is a notable increase in gene-derived piRNAs toward the end of the cycle. Our results suggest regulatory roles of piRNAs in modulating both TEs and mRNAs during BOF maturation, likely related to changes in the follicular cell program.

Ancient human dental calculus is one of the richest archives of archaeological biomolecular information, providing direct evidence of diet, oral health, and the oral microbiome. Proteomic analyses of this biological matrix have so far focused mainly on oral microbes and dietary proteins, with milk proteins such as beta-lactoglobulin (BLG) providing the largest corpus of proteomic evidence. Despite the close relation between the various stages of dental calculus formation and mineralization with the dental enamel surface, proteins from the dental enamel matrix have not previously been reported outside of dental enamel tissue. Here we reanalysed 498 ancient dental calculus proteomes from 14 published studies (n=434 individuals) reporting the presence of BLG, spanning from the Neolithic to the Victorian Era and applying different protein extraction protocols (FASP, GASP, SP3 and in-solution digestion). Dental enamel matrix proteins were identified in ten studies (n=37 individuals), with amelogenin being the most frequently detected. Enamel peptides occurred more often in studies that applied SP3, although amelogenin was successfully identified through both SP3 and FASP. Structural proteins, including enamelin, ameloblastin, and MMP20, were also identified. The detection of AMELX and AMELY peptide sequences provided new insights into cases where the sex was previously undetermined. These findings establish dental enamel proteins as a new category of biomolecules detected in dental calculus, broadening its application beyond diet and microbiome studies to possible sex estimation. HighlightsO_LIDental calculus entraps oral microbes along with endogenous and exogenous particles during formation and mineralization C_LIO_LIWe conduct reanalysis of 14 published ancient dental calculus studies (n = 434 individuals) spanning the Neolithic to Victorian Era C_LIO_LIDental enamel proteins AMELX, AMELY, AMBN, COL17A1, ENAM and MMP20 are identified in ancient human dental calculus C_LIO_LIAmelogenin was the most frequently detected enamel protein C_LIO_LIWe expand dental calculus palaeoproteomics beyond diet and oral microbiome to potentially include sex estimation C_LI

Medusozoan cnidarians (e.g., jellyfish) metamorphose from a benthic juvenile polyp into a pelagic adult medusa, providing a well-known example of a clade that uses tissue remodeling to create distinct juvenile and adult body plans. Staurozoans (i.e., stalked jellyfish) are an atypical lineage of medusozoans that have lost their medusa stage; thus, their juvenile and adult body plans look remarkably alike. Their limited metamorphosis is characterized by the regression of primary (juvenile) tentacles and the development of secondary (adult) tentacles. In some staurozoan lineages, metamorphosis also involves development of novel adhesive structures (anchors), which are built on top of the regressing primary tentacles. Understanding how cells are partitioned from making juvenile tissues to making adult tissues is important for understanding how animals can make adult structures in the absence of complete metamorphosis. We compared the abundance and distribution of proliferative cells in tissues undergoing regression (primary tentacles) and development (secondary tentacles and anchors) during the juvenile to adult transition in the San Juan Island stalked jellyfish, Haliclystus sanjuanensis. We show that proliferative cells are lost in regressing primary tentacles but are gained in anchors, consistent with a shift in investment from juvenile to adult tissue. Prior to regression, primary and secondary tentacles show similar patterns in their proliferative cell distribution and in the identity of their cnidocytes (stinging cells), indicating that adult tentacles are made by re-deploying a juvenile tentacle program. Finally, we demonstrate that unlike secondary tentacles, primary tentacles cannot regenerate, illustrating that the temporary investment in this tissue is tied to their loss of proliferative cells. Thus, we propose that continued investment in a population of proliferating cells is an important mechanism for segregating temporary tissues (primary tentacles) from long-term tissues (secondary tentacles). These observations of cell dynamics in H. sanjuanensis suggest that temporary investment into juvenile structures may be used to pattern novel adult tissues, providing an important mechanism for diversifying adult body plans.

Life history traits are often correlated, creating trade-offs that may impede the response to natural selection and be responsible for the evolution of senescence. These trade-offs may arise through pleiotropic effects, which can affect the response to selection in ways that resemble intra-locus sexual antagonism. Despite these hypothesized relationships, we lack clear connections between pleiotropy, sexual antagonism, and the evolution of life histories. Empirical tests for inter-sexual differences in life-history traits, including sex-specific aging, can be used to evaluate hypotheses about how pleiotropy and sexual conflict affect evolutionary trade-offs. To those ends, we measured lifespan, development time, and body size in Drosophila pseudoobscura males and females, each of which carried one of six third chromosome inversion genotypes. Temperature affected lifespan and development more than any other factor; higher temperatures increased mortality rate, decreased lifespan, and accelerated development. However, we also observed sex differences in mortality rates and development times that depended on genotype and temperature. Notably, temperature elevated the initial mortality rate across all flies, yet increasing temperatures reduced the rate of aging in some genotype-sex combinations. Similarly, direct effects of genotype on mortality rate and development time depended greatly on sex and temperature, but there was no genotype effect on body size. Despite these context-dependent genotype effects on life history traits, we failed to identify any correlations that would serve as clear evidence for sexual conflict or trade-offs. Our results therefore suggest that either historical conflicts have been resolved or any conflicts that may exist do not result in the correlations predicted by existing models.

The evolution of visual systems has compelled numerous investigations of developmental processes underlying eye patterning across Bilateria. It is well-established that homologs of the transcription factor Pax6 play a highly conserved role in eye fate specification and are at the top of the retinal determination gene network (RDGN) hierarchy. In insects, the two Pax6 homologs eyeless (ey) and twin of eyeless (toy) are required for the development of the two visual systems broadly found within the phylum (i.e., median and lateral eyes). Curiously, Pax6 homologs do not appear to maintain this function in well-studied chelicerate models, with emphasis on spiders, a lineage of arachnids with great diversity of eye form and acuity. It was recently proposed that the gene Pax2 (shaven; sv) may have subsumed the role of eye fate specification in chelicerates, a hypothesis predicated upon the observation that one of two spider Pax2 copies is strongly expressed in the developing lateral eyes during embryogenesis. However, no functional data are available for any Pax homologs across Chelicerata. We examined the incidence of Pax family genes across Chelicerata, as well as interrogated the expression and function of Pax2 and Pax6 homologs in the daddy-longlegs Phalangium opilio, an arachnid recently discovered to bear a highly plesiomorphic arrangement of visual systems. Here, we show that ey and toy are expressed early in the developing head lobes of P. opilio, whereas sv is not expressed until well after stages when downstream RDGN members (eyes absent and sine oculis) are already activated. Gene silencing of ey, toy, and sv individually had no discernible effect on eye development. By contrast, double knockdown of ey and toy resulted in an array of median eye defects, spanning loss of some cells of the eye to total loss of the median eyes. Gene expression assays also showed that depletion of the two Pax6 copies resulted in failure of the vestigial median and vestigial lateral eyes. These data are consistent with a conserved role for Pax6 homologs in patterning both visual systems and all three eye pairs in the daddy-longlegs. Our results comprise the first functional data for Pax6 genes in any chelicerate and suggest that heterochronic shifts in expression, rather than changes in function, underlie the atypical dynamics of Pax genes in derived arachnid groups such as spiders.

Animals sense and integrate complex external cues to make developmental decisions that help them better survive and adapt to their natural habitats. Under environmental adversity, nematodes can enter an alternative developmental pathway to form a diapautic and stress-resistant stage, termed the dauer larvae. While dauer formation has been well characterized in Caenorhabditis elegans, how environmental factors influence analogous stages in other nematode species remains largely unexplored. This study examines how symbiotic bacteria, temperature, and pheromones affect the formation of the infective juvenile (IJ), a dauer-like stage, of the insect-parasitic nematode Steinernema hermaphroditum. In contrast to C. elegans, where dauer entry is promoted by heat, IJ development in S. hermaphroditum development is enhanced by reduced temperature. Moreover, the presence and absence of live symbiotic bacterium Xenorhabdus griffiniae functions as an ON-and-OFF switch that regulates the host IJ formation. Crude pheromone extracts from S. hermaphroditum liquid culture do not robustly induce IJ formation in a dose-responsive manner, unlike the potent pheromone-driven dauer entry observed in C. elegans. Nutrient-rich liver-kidney media that mimics host insect environment showed IJ entry induction in a pheromone-dependent manner. These data suggest that external cues, such as temperature, microbial diet, and pheromone, are perceived differently by S. hermaphroditum in comparison to that of C. elegans, reflecting species-specific adaptations to distinct ecological niches and life history strategies.

Spatiotemporal organisation of biological molecules is a key driver of cellular processes, including many post-transcriptional epigenetic processes. The germline-specific germ granules are biomolecular condensates that act as hubs for mRNA and small RNA processing and are core regulators of germline gene expression programming. Germ granules have been studied extensively in C. elegans, and recent developments have led to many subdivisions of the germ granule into specialised compartments. Rapid advancements in microscopy and protein-protein interaction (PPI) screening techniques have produced a large amount of data towards characterising the localisation of proteins to specific granules. However, common methods used to probe PPIs are limited in their ability to robustly detect valid interactions, especially the multivalent and sometimes transient ones observed in granule environments. Here we perform a meta-analysis of granule protein interaction screens. While these experiments generally enrich for proteins matching the profile of granule-associated proteins, we find that when considering screens individually, reproducibility is surprisingly low, highlighting not only the variability inherent in these methods but also the dynamic nature of the PPI networks present in granules. We developed an algorithm to provide a measure of each proteins association with specific granules across various experiments. By further clustering and investigation of the resulting score matrix, we demonstrate the power of this holistic approach to provide deeper insights into germ granule organisation and highlight novel can provide a resource to better inform future investigations into granules and their constituent proteins.

Parental rejection of apparently healthy newborns is widely classified as a behavioural abnormality in captive primate colonies, yet its biological significance remains unclear. In owl monkeys (Aotus nancymaae), parental rejection, defined here as cessation of nursing leading to rescue nursery rearing, is typically lethal for offspring and is transmitted across generations despite reducing offspring survival. Here, we tested whether parental rejection is associated with lifespan and reproductive differences in parents and their surviving offspring. We analysed long-term demographic records from a captive colony of 962 individuals and compared survival and reproductive outcomes between rejector and non-rejector parents using survival analyses and regression-based models. Parents who rejected offspring lived significantly longer than non-rejectors, with an average lifespan advantage of approximately 4-4.5 years in both males and females. This survival difference was concentrated during the prime reproductive period (6-20 years). Well-reared offspring of rejector parents also lived longer than offspring of non-rejectors, with a mean lifespan difference of 1.26 years. Rejector parents produced more offspring overall, but this difference was explained by extended lifespan rather than higher reproductive output per year. Analyses stratified by rejection timing showed no longevity advantage in first-birth rejectors, whereas parents rejecting later-born offspring exhibited longer survival. Together, these findings show that parental rejection is associated with longer lifespan in parents and in their well-reared offspring under captive conditions. These patterns are consistent with altered allocation of parental investment under energetic or environmental stress.

Conspicuous coloration in animals is generally thought to evolve and be maintained through inter- or intraspecific interactions such as mate choice, but this might not always be the case. The sight-line hypothesis proposes that conspicuous light-dark contrast in front of the eyes (hereafter, eyeline) evolves and is maintained due to viability selection, enhancing an individual visual acuity and thus evolutionarily associated with a particular foraging behavior that requires accurate aiming. However, empirical evidence that supports the sight-line hypothesis is virtually absent, with no studies demonstrating the key prediction that the direction of eyelines matters. Here, I tested the sight-line hypothesis using macroevolutionary analyses in terns and allies, which are a suitable study system, because they have variation in facial color patterns, including presence/absence and, if any, various angles of eyelines. They also have a large variation in foraging behavior, including picking, plunge diving, and skimming. As predicted by the sight-line hypothesis, tern lineages that require accurate aiming at foraging (e.g., plunge diving) are more likely to have eyelines. In addition, the evolutionary transition to the state with eyelines and these foraging behaviors was more likely to occur than the reverse transition. Furthermore, as expected by the fact that the direction of travel is upwardly deviated from the direction of the bills during skimming, the eyeline angle from bills was evolutionarily positively associated with the occurrence of skimming behavior. To my knowledge, the current study is the first to demonstrate that the direction of the eyeline matters, thereby strongly supporting the sight-line hypothesis.

Phenotypic diversity is often thought to arise from the evolutionary modification of developmental processes. However, developmental processes are tightly coupled in space and time, with each process beginning from conditions set by the one before it. While we know from dynamical systems theory that initial conditions can significantly affect a systems out-come, their importance as a source of phenotypic evolvability has been largely overlooked. Here we show for the first time, that phenotypic evolution can proceed through changes in developmental initial conditions while the underlying developmental process remains conserved. Somitogenesis is the process by which vertebral precursors, known as somites, are periodically patterned in the pre-somitic mesoderm (PSM). Somitic count (total number of somites) is thought to diversify through the evolution of components of somitogenesis such as the tempo of the segmentation clock or the mechanisms driving axial morphogenesis. Using two closely related species of Lake Malawi cichlid fishes that differ in vertebral counts, we show that somite count evolution has happened without changes to somitogenesis itself, but instead, by altering the size of the PSM at the onset of this process. This work will expand what we consider developmental drivers of phenotypic evolution and highlight the importance of comparative studies to understand the diversification of phenotypes.

We present the first radiographic ageing framework for common dolphins (Delphinus delphis), based on ossification and epiphyseal fusion patterns in the pectoral flipper, demonstrating higher reliability for chronological age estimation than currently available epigenetic approaches for this species. Using individuals of known dental age, we calibrated two modelling approaches to predict dental age from radiographic bone scores: 1) a univariate polynomial regression using a total bone score (sum of 16 scores across all assessed flipper bones), and 2) a multivariate canonical analysis of principal coordinates (CAP) incorporating 16 individual bone-score variables. Both approaches successfully predicted dental age from skeletal ossification patterns. For an age range of 0 to 24 years, polynomial regression demonstrated high predictive accuracy with median absolute errors (MAEs) of 1.25 years in females (Spearmans {rho} = 0.93, R{superscript 2} = 0.90) and 1.08 years in males ({rho} = 0.95, R{superscript 2} = 0.86). The CAP model yielded MAEs of 1.35 years in females ({rho} = 0.90, R{superscript 2} = 0.85) and 1.80 years in males ({rho} = 0.94, R{superscript 2} = 0.84). Notably, both radiographic bone ageing models achieved equal or lower median absolute errors and higher coefficients of determination than a recently developed epigenetic clock for common dolphins derived from the same population (MAE = 1.80, Pearsons correlation (r) = 0.91, R{superscript 2} = 0.82). When applying the bone ageing models to individuals of unknown dental age, both models produced age estimates consistent with expected life-history stages (foetus, neonate, juvenile, subadult, adult), although accuracy declined in dolphins above 20 years, likely as a consequence of subtle age-related variation in skeletal changes in this species. Radiographic ageing provides an accurate non-invasive tool for demographic assessment to support conservation management of common dolphins.

Our bodies have evolved to maintain homeostasis through regulatory systems that continuously adapt to keep physiological processes within a normal range. From this perspective, complex chronic disease can be understood as a breakdown of compensatory mechanisms, resulting in loss of homeostasis. Here we propose that adaptive receptor expression dynamics may serve as one such compensatory mechanism, increasing receptor surface expression when external ligand is insufficient, and clearing it when signaling is excessive. To explore this, we adapt a previously published agent-based model and simulate it under a range of scenarios. We find that the system of adaptive receptor expression is robust to oscillatory perturbations but not to chronic stress. We propose that receptor turnover dynamics may be better understood as an adaptive, environmentally responsive process rather than a fixed biological property, and that in some cases, disease manifests only after compensatory mechanisms have been pushed past their limits. We conclude with a discussion of implications for understanding complex chronic diseases, for thinking about epigenetic and mutational change as escalating layers of adaptation, and for how we model receptor dynamics in the context of receptor-mediated drug activity.

Transcriptomic analyses are widely used to elucidate the molecular mechanisms driving gametogenesis and reproduction in fish, yet their accuracy depends heavily on appropriate normalization of gene expression data. Conventional approaches that rely on single or multiple reference genes are problematic during teleost oogenesis, as profound structural and physiological remodeling of the ovary challenges the assumption that commonly used reference transcripts remain stable. In this study, we assessed by qPCR the transcriptional variability of four widely used reference genes (actb, ef-1, gapdh, and 18S rRNA) throughout the oogenic cycle of the thicklip grey mullet (Chelon labrosus), using geNorm and NormFinder analyses, and we additionally evaluated total cDNA concentration as an alternative normalization factor. To examine the performance and interpretive consequences of each normalization strategy, we compared expression patterns of key steroidogenic genes (star, cyp19a1a, and cyp11b) normalized by individual reference genes, combinations of reference genes, or total cDNA concentration. All evaluated reference genes displayed notable transcriptional variability across oogenesis, confirming their limited suitability as sole internal controls. In contrast, normalization approaches integrating multiple reference genes and/or total cDNA concentration consistently provided greater stability and more reliable biological interpretation. These results support a refined and more robust normalization framework for transcriptional analyses in fish ovaries, particularly during stages of extensive tissue remodeling. Our findings demonstrate cDNA-based normalization is straightforward, rapid, and easy to implement across laboratories, providing a practical alternative for achieving accurate, reproducible transcript quantification in fish ovary studies.

Ehmt2 is a key H3K9 methyltransferase that regulates genome silencing and structural integrity during animal development. In addition to this canonical function, Ehmt2 has also been implicated in neural tissues mediating both direct and indirect transcriptional activation, and exon splicing, to facilitate proper neural cell differentiation and survival. Several germline loss-of-function animal models have been developed showing both conserved and divergent phenotypes that range from embryonic lethality to behavioural deficits in adult, fertile animals. Here, we generated the first maternal-zygotic ehmt2 loss of function mutant in zebrafish using CRISPR-Cas9 mutagenesis. An assessment of the pattern of H3K9 methylation in mutant embryos by ChIP-seq indicates that there are aberrant levels of this repressive mark, including reduction in discrete 5 non-coding regions of genes, but with no significant change in the overall pattern distribution of these marks across the genome. Global transcriptome and morphological analyses demonstrated that mutant embryos displayed greater variation in the timing of developmental progression that is, on average, slower compared to controls. Despite this, mutant embryos ultimately survive and are fertile. Through examination of progenitor cell dynamics and gene expression profiles, we found that the delay in embryonic development was associated with longer rates of S-M phases of the progenitor cell cycle in mutants leading to deficits in tissue growth. Finally, our data suggest a robust network of epigenetic regulators can potentially compensate for Ehmt2 loss of function and permit embryonic development and survival in ehmt2 mutant zebrafish. Our work establishes a zebrafish ehmt2 loss of function model that will facilitate examination of the complex and varied roles of Ehmt2 in vertebrate development.

Sex-responsive trait development generates much of the phenotypic variation found in natural populations and diversifies rapidly among closely-related taxa. Furthermore, rather than exhibiting equal sexual dimorphism across all traits, organisms are mosaics of tissues that vary in their degree of dimorphism. Yet, how these mosaic patterns are generated remains largely an open question, as sexually dimorphic traits have typically been studied individually in select model systems. In this study, we compare gene regulatory landscapes across five traits that differ in the degree of morphological sexual dimorphism in the bull-headed dung beetle Onthophagus taurus by assaying tissue-specific gene expression and chromatin accessibility at the onset of pupal development when future adult form is specified. We identify a modest number of pleiotropic regulators associated with sex differences across traits, yet uncover a high degree of sex- and trait-specificity in chromatin architecture within developing tissues. We then confirm the role of the sex determination factor doublesex in the regulation of sex differences through expression of sex-specific isoforms, and uncover trait- and sex-specific sets of Doublesex binding sites likely underpinning context specific sexual dimorphisms. Further, we identify and functionally validate the transcription factor ventral veinless as a regulator of sexually dimorphic development. Our findings suggest that in contrast to doublesex, ventral veinless does not exhibit sex-biased expression, yet exerts its sex-specific regulation via sets of differentially accessible binding sites. This work furthers our understanding of the molecular mechanisms instructing the development of sex differences and provides novel insights illustrating how transcriptional activity and chromatin remodeling interact to generate sexual dimorphism in a trait-specific manner. More generally, our work contributes to a growing body of knowledge on how development integrates cues such as sex determination to enable highly similar genomes to yield diverse phenotypic outcomes.